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Abstract 

Bohmian trajectories on the toroidal surface are determined from eigenfunctions of 
the Schrodinger equation. An expression for the monodromy matrix M(t) on a curved 
surface is developed and eigenvalues of M{t) on calculated. Lyapunov exponents 
for trajectories on are found for some trajectories to be of order unity. 
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The equation of motion for calculating Bohmian trajectories follows from inserting \E'(r, t) = 
R(r,t)e^^^^'*^ into the Schrodinger equation and separating the ensuing expression into real 
and imaginary parts [1]. The result is {m = h = 1) 

'i--V(V^Q) (1) 
with V = VS and Q the quantum potential 



Q = - 



2 R 



Given a distribution of initial positions as determined from ?/'*?/' the standard results of 
quantum mechanics are recovered. 

This letter concerns the influence of curvature on Bohmian trajectories. Wu and Sprung 
[2] have pointed to the importance of a trajectory's initial position as critical to its evolution 
in time. If the trajectory lies on a curved surface, local curvature near the trajectory's initial 
position may be expected to play a role in its character. Here wave functions on a curved 
surface are employed to generate trajectories from S [2,3,4]. 

was chosen as the surface with which to investigate curvature effects on Bohmian 
trajectories for several reasons: First, the torus has non-trivial mean and Gaussian curva- 
tures[5]. Second, good approximate wave functions for a particle on are available [6]. 
Finally, the rectangular strip R^xR^ with periodic boundary conditions can be used as a flat 
torus analog for comparison between trajectories on it and those on T^. For convenience the 
strip will be referred to as F^. 
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A toroidal surface with major radius R and minor radius a may be characterized by the 
Monge form 

r(^, (f)) — {R + a cos9)ep + a sinO (3) 

for which 

ds^ = a^de^ + (R + a cosefdcf)'^. (4) 

The Hamiltonian is taken as H — — |V^. Defining a — (3 — 2Ea?, and making the 
standard x(0) = e'"*"^ ansatz for the azimuthal eigenfunction gives the Schrodinger equation 

5^-0 a sin 6 dip m^a^ i r i r\ fK\ 

W~ [l + acose]de~ [1 + a cose]^"^^ " ^ ' 

Eq. (5) was solved in [6] by a Fourier method. The explicit forms of the surface toroidal 
wave functions (STWs) used here are given in table I for R — l,a — 1/2. The eigenfunctions 

on analogous to those in tabic 1 arc ( |e*™''^. 

As evidenced in table 1, six STW states were chosen to comprise three positive parity 
states and three negative parity states. The [nm] values selected for the negative parity 
states arc the same as those selected for the positive parity states. This [nm] matching 
was motivated by a desire to construct (as closely as possible) functions for comparison to 
gine functions on F^. Additionally, without superpositions of the form -0+^ ± some 
interesting motion about the minor radius of would not be manifest. 

The phase 5" for and is generated with a six state superposition 

^{e, 0, t) = E Cnm'^nmie, 0)6"^^-*. (6) 

Trajectories were determined from 



nm 



S = tan-^ 

and 



(7) 



Eqs. (7) and (8) yield many classes of surface trajectories. Figures 1 and 2 give two 
state results for to the left of each figure and for to the right. Figure 3 shows a 
quantized trajectory structure that emerges from a combination of negative parity states 
on T^. Figure 4 is a (0,0) plot for a path on showing rapid variation in phase space. 
Figure 5 is an example of a {6, 6) plot for a superposition with a dominant mode and small 
admixtures of two other states. A small change on the order of a few parts in 10^^ in 
causes ample modification to the path. F^ phase space plots have not been shown because 
they demonstrate (at least for the cases above) very little structure. 
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A measure of the divergence of two initially nearby trajectories is associated with the 
eigenvalues of the monodromy matrix M{t) [7,8], used to good effect by Frisk [3] for the 
study of Bohmian paths on a two-dimensional rectangular geometry. M{t) gives the time 
evolution of the separation vector t^x via 

(5x(i) =M(t)5x(0). (9) 

Since 

j^[S^{t)]^V[5S], (10) 
the monodromy matrix can be shown to obey the differential equation 

JM. (11) 



dt 

For R^xR^ with or without periodic boundary conditions the J matrix is [3] 

j^fS S^y\ (12) 

On a two dimensional surface characterized by a locally orthogonal metric 

ds^ = Quudu^ + gvvdv^ 
Eq. (12) must be extended to (written in a cumbersome but illustrative form) 

On T2 Eq. (13) yields 

( See ^e± \ 
aG G2 / 

with G = l-\- a cos 9. 

It is worth noting that for curved surfaces 5 and V do not commute. The ordering chosen 
in Eq. (10) insures all quantities he on the local tangent space [9,10]. 

The differential equations implicit in Eq. (11) were found to be time consuming and slow 
to converge to suitable accuracy with standard methods on T^. Further, it was not clear 
that the results were accurate to any order. However, the foremost goal here was obtaining 
comparisons between the eigenvalues of M{t) on versus those on for several values 
of Oq. With this is mind, a Mathematica code was set to solve the differential equations 
to a lower accuracy for the relatively short time t — 10. The advantages of adopting this 
procedure were a) confidence that our numbers were accurate to at least four significant 
digits and b) each point took only at most two minutes to acquire on an 800 MHz rated 
home PC [11]. 
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Rather than work directly with the eigenvalues (5i of M{t), a Lyapunov exponent Aj is 
defined through the relation 

Xi{t)^jhiPi{t). (15) 

A more convenient definition which serves to dispense of an unknown overall norm is one 
employed in quantum Monte Carlo calculations [12] and adopted here, 

A. = ^^l^#4- (16) 

The larger of the two values of Aj = A was taken. 

Tables II and III give values of A as a function of 9q for two sets of c„„j values. For both 
sets A is generally an order of magnitude larger on than on F"^ (save for one anomalous 
point). The F"^ values show a symmetry about 6 = it while the values do not. The 
results certainly show dependence on 9o but do not show any simple relation to the Gaussian 
or mean curvatures K and if on T^. This is not unexpected. The Hamiltonian of Eq.(5) 
does not incorporate H and K explicitly but rather factors related to them in a complicated 
manner. A trajectory on generated from a representative point in table II is given in 
figure 6, and one from table III is shown in figure 7. Again, the corresponding F^ plots do 
not possess enough structure to warrant their inclusion. 

In this letter linear combinations of wave functions were employed to generate surface 
and phase space plots of Bohmian trajectories on and F^. The plots and tables illustrate 
that curvature can alter trajectory structure. This is a direct manifestation of the form of the 
wave functions on a curved surface. Each exponential function in a superposition of states 
on possesses its own time dependent energy phase while each STW is a superposition of 
several functions attached to one phase. 

Cnm combinations making Lyapunov exponents on appreciably smaller than those on 
F^ over many values of Oq have not been found. Positive A of order unity is usually taken 
to indicate chaotic behavior, but while some trajectories are certainly complex, and A an 
adequate measure of the distortion of local tangents, we consider it premature to state the 
trajectories are chaotic [13,14]. Nevertheless, it has been shown that a few low-lying states 
can yield complex phase space behavior on the the torus. Because the torus is a simple 
compact surface, its curvature likely causes less modification to trajectory structure than 
surfaces with rapidly varying regions of H and K. Those surfaces are certain to induce 
greater complexity in the trajectories. 
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Figure captions 



Fig. 1. Trajectories for the superposition ^(^,0, t) = ^§^'32 + ^y|^32 with t — 30 and 
^0 = 0. The trajectory appears to the left of the figure. The trajectory projected onto 
the torus is to the right. 

Fig. 2. Trajectories for *(6', 0, t) = ^32 + «V^*2i with t = 24 and 60 = I.OStt. 

Fig. 3. Trajectories for ^{6, (p, t) = ^"^ti + (\/W*2i - ^•\/M)^2i with t = 50 and 60 = 0. 

Fig. 4. ^) plot for 0, = ^/l'^2l + ^yi^32 - yi^32 with t = 30 and 0o = 1.257r. 

Fig. 5. (6,6) plots for = y^^ro + V^*2i + /i^32 with t = 30 and 

00 — 1.4247r shown on the left and for Oq — 1.4297r shown on the right. 

Fig. 6. {6, 6) plot for ^(^, 0, t) = ^1^+ + ^^32 with t = 36 and ^0 = 0. 
Fig. 7. 0) plot for (/), t) = + «v7*ro with t = 36 and ^0 = 0. 
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Table I: Surface torodial wave functions and eigenvalues for R — 1, a — 1/2. Coefficients not 
listed are at least an order of magnitude smaller than those given. 



*±^;i?=l,a = l/2 




^ro = 


-.2176 + .4352 cos^ - .0714 cos2^ + .0118 cosS^ 


1.2223 


^2+1 = 


-.0733 + .2419 cos^ - .8393 cos2^ + .0541 cos3^ 


4.4767 


*3+2 = 


-.0420 + .1240 cos9 - .2772 cos2^ + .8240 cos3^ 


10.6657 


*ro = 


+.8118 sin^ - .0739 sin2^ 


0.9767 


*2"l = 


-.1799 sin^ - .8367 sin2^ 


4.4106 


*3"2 = 


-.0808 sine + .2568 sin2e - .8257 sin3e 


10.6151 



Table II: Lyapunov exponents Ag measured ait — 9, Aio measured at t = 10 and A as defined 
by Eq. (16) for the superposition ^1^32 + \l\'^^2- The upper half of the table gives results 
for and the lower half are results for F^. 
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TT 
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TT 
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27r 
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TT 
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47r 
3 


37r 
2 


b-K 
3 


IItt 
6 


Ag 


2.12 


4.23 


2.53 


1.91 


3.36 


2.29 


2.75 


2.55 


3.11 


1.72 


1.23 


2.70 


Aio 


4.10 


5.08 


3.67 


2.65 


4.17 


2.69 


2.57 


2.85 


5.73 


2.39 


2.59 


3.27 


A 


21.9 


12.7 


13.9 


9.35 


11.5 


6.35 


.96 


5.53 


29.3 


8.44 


14.8 


8.40 


Ag 


.030 


.036 


.034 


.033 


.036 


.032 


.030 


.036 


.034 


.033 


.033 


.036 


Aio 


.047 


.056 


.050 


.046 


.049 


.048 


.047 


.056 


.050 


.046 


.046 


.056 


A 


.185 


.234 


.179 


.161 


.155 


.195 


.186 


.233 


.179 


.160 


.155 


.235 



Table III: Lyapunov exponents Ag measured at i = 9, Aio measured at i = 10 and A as 
defined by Eq. (16) for the superposition ^J\^to + \l\^io- The upper half of the table gives 
results for and the lower half are results for F^. 
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3 
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Ag 


2.46 


1.69 


1.68 


1.63 


1.67 


1.86 


2.01 


3.59 


3.64 


3.22 


3.13 


6.87 


Aio 


2.55 


1.65 


1.63 


1.58 


1.61 


1.74 


1.92 


3.28 


3.97 


3.99 


3.73 


6.90 


A 


3.39 


1.23 


1.02 


1.19 


1.07 


.067 


1.18 


.447 


7.00 


10.90 


9.14 


7.16 


Ag 


.086 


.684 


-.003 


.022 


.002 


.009 


.086 


.684 


-.003 


.022 


.002 


.009 


Aio 


.122 


.336 


.450 


.124 


.015 


.032 


.122 


.336 


.450 


.124 


.015 


.032 


A 


.415 


-2.45 


4.52 


.942 


.133 


.239 


.415 


-2.45 


4.52 


.942 


.133 


.239 
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Figure 1: 




Figure 3: 
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Figure 6: 




Figure 7: 
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